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Abstract 

We apply the method of zeta functions, together with the n* -prescription for the tem- 
poral gauge, to evaluate the thermodynamic pressure in QCD at finite temperature 
T. Working in the imaginary-time formalism and employing a special version of the 
unified-gauge prescription, we show that the pure-gauge contribution to the pressure 
at two loops is given by P| auge = -(g 2 /lH)N c N g T A , where N c and N g denote the 
number of colours and gluons, respectively. This result agrees with the value in the 
Feynman gauge. 
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1 Introduction 



The temporal gauge is a physical, i.e. ghost-free gauge which belongs to the class of axial- 
type gauges characterized by the constraint 

n^(x) = 0, = 0,1,2,3,..., (1) 

where A^(x) is a massless Yang-Mills field, and a = 1,2, . . . , N 2 — 1, for SU(N). These 
axial- type gauges also include the pure axial gauge [n 2 < 0), the planar gauge (n 2 < 0), 
and the light-cone gauge (n 2 = 0). In the temporal gauge, the fixed vector n M is time-like: 
n 2 = Uq — n 2 > 0. Today, the temporal gauge is as useful as it was 65 years ago in the 
quantization of the Maxwell- Dirac field by Weyl in 1929 U and Heisenberg and Pauli in 
1930 @. Only the degree of complexity has increased: today's temporal gauge is applied 
in such sophisticated cLl CclS clS the vacuum tunneling by instantons ||, Nicolai maps j|, [|], 
one-loop thermodynamic potentials , and especially in quark-gluon plasma studies JJJ at 
finite temperature T. Of course, there are good reasons for this popularity. 

Due to the space-time asymmetry at finite temperature in the imaginary-time for- 
malism, the thermodynamic equations are no longer symmetric in 4-space, but are only 
invariant under spatial rotations. Accordingly, since the temporal gauge-fixing condition 
noA^x) = leaves the spatial components of A a ^(x) unconstrained, the temporal gauge 
may indeed be regarded as the "appropriate gauge" for this particular boundary-value 
problem. 

During the last ten years, finite-temperature calculations in the temporal gauge have 
been performed both in the real-time formalism and in the imaginary-time formalism. 
Kapusta and Kajantie f7], for instance, employed the imaginary-time formalism to evaluate 
the leading terms in the 1/T-expansion for the two-point function (static limit). The finite- 
temperature case was also examined by Landshoff and James |§, as well as by Brandt, 
Frenkel and Taylor |J. Working in the Feynman gauge, Brandt, Frenkel and Taylor utilized 
the method of ^-functions to derive the complete 1/T-expansion for the one- loop Yang-Mills 
self-energy. 

In 1994, the present authors developed a general procedure for doing perturbative 
calculations in the temporal gauge at both zero (T = 0) and finite temperature (T ^ 0). 
The procedure hinged on a special version of the n*-prescription, originally designed for 
the spurious poles of the temporal-gauge propagator (cf. Eq. (^)). No problems were 
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encountered for zero temperature, but for T^Oa different approach was needed. In order 
to simplify the computation of Matsubara frequency sums, it became necessary to replace 
the traditional contour method by the method of zeta functions. The latter method is 
generally known to facilitate the computation of high-temperature expansions to any order 
in 1/T. 

Application of the ^-function technique, together with the n*-prescription for the tem- 
poral gauge, enabled us to evaluate the complete 1/T-expansion for the self-energy Yl^k^ = 
0, k) |TTJ, and to derive the Debye chromo-electric screening length m e j in the infrared limit, 



The purpose of the present article is to apply this successful procedure [[U| to another 
important quantity, namely the thermodynamic pressure. Working in the imaginary-time 
formalism, we shall evaluate the thermodynamic pressure at finite temperature to one and 
two loops. 



2 Basic Tools 



The temporal-gauge propagator for Yang-Mills theory at finite temperature reads as follows 



G±(P) 



-i5 



ab 



(2tt) 2uj (p 2 + ie) 



_ p»n u + p u n^ n 2 p il p h 



n 2 > 0, e > 0, 



p ■ n (p • n) 2 

where 2u denotes the dimensionality of complex space-time; in four dimensions, g^ v 
diag(+l, -1, -1, -1), /j,, v — 0, 1, 2, 3. For = (1, 0, 0, 0), the propagator reduces to 

-i5 ab 



(2) 



(2tt) 2w (p 2 + ze) 



- 6a + 



PiPj 



(Pb) 



i,3 = 1,2,3, 



GS(P) = G%{p) = G%{p) = 



tab 



ab I 



(3) 



The spurious double pole at p = in Eq. (|3|) may be treated by a special version of 
the unified-gauge prescription for axial- type gauges JTB], [TJJ, namely 



1 

Po 



temp 
temp 



Po 



lim 

Po + ie ' 



e > 0: 



(Po 



lim ( - 

\q 



Po 



e > 0. 



(4) 



"Po + it 

Note that prescription (§) is causal and well defined for all p = 2ninT, n = 0, 1, 2, . . ., T 
being the temperature. In particular, no ambiguities arise for the case n — 0. 
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In addition to the formulas (Q) and (|j), loop calculations at finite temperature also 
require the replacement 

d 4 p . ^ n l^ f rf 3 p 



(2vr) 4 /3 n f^ (2tt)3- 
Here (3 = l/(kT), and = iw„, with u; n = 27rnT; it is customary to equate Boltzmann's 
constant k to unity. Finally, we note that the components of n M = (no, ni, ri2, n$) have the 



following structure in the temporal gauge 



so that 



= (no, n ± , -i|nj.|), n ^ 0, n ± = (m, n 2 ); 
n* = (no,njL,+i|n-i|)> ( 6 ) 



lim n M = (1,0,0,0), n = l. (7) 

|nj_|-+0 



3 The Q.C.D. Pressure to One Loop 

Since the pressure P is defined in terms of the partition function Z, 

P = T^— — , T = temperature, V = volume, (8) 
oV 

we first need to calculate Z for pure gauge theory [f[5 



Zgauge = / [dA„] det(n ■ d)e s ; (9) 



denotes a product over all field configurations, with the action S given by 



o 



S = I dr I d 3 x 



2_ 



4 2a 



= \[dr\SxA^(dH, v -d,d v -^)A\ (10) 

We shall keep n M general for now (cf. Eqs. (M), taking the limit |nj_| — > at a later time. 
Notice that periodicity of the boundary condition causes the boundary term in Eq. ( P"U|) to 
disappear. Introducing the Grassmann variables C and C (C is the hermitian conjugate of 
C), we may write the determinant det(n • d) as follows: 

det(n -d)= J [dC)[dC] exp (J dr J d 3 xCd-nC). (11) 
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It is convenient now to apply a discrete Fourier transformation in Euclidean space (the 
continuous limit will be taken later in the calculation): 

^,r) = -|rEE e,(p,X+U " r> ^(p). "n = (2n + 1)ttT, (12) 

V 2 „ p 

where C n (p) is the Fourier transform of C(x, r), p and n label the momentum and energy- 
quantum numbers, respectively, and where the volume factor V is necessary for proper 
normalization of the transformation. The frequency u n is odd here, because C is a fermion 
field. For the gauge field A M (x, r), we have: 

A^^r) = M^T, el{P '^ nT)A n(p), uj n = 2nnT. (13) 

' " n p 

The extra factor of a//? is needed to keep the action dimensionless. Thus 
Zgauge = J [dA"] [dC] [dC] exp { J dr J d 3 xCn ■ dC) 

exp{ f dr [ d 3 /EE e^« + ^')V(P + P'>M£(pKV) 

J ° J ZK n,Pn',p' 

(-P 2 ^+P^-^)}, (14) 
= f [tL4"] [dC] [dC] exp I f dr f d 3 x^- ]T ]T e -*( P -x+^r) 

J { J ° J V n,p n',p> 

a n (p)(ip0e^'- x+ ^)(7 n Kp0}exp|^(WE^(p)^-n(-p) 
(-P 2 ^+P^-^)}, (15) 

= J [dq[dC']exp|^5:d„(p)(^p)e^ x+ ^)C' n (p)| 

(16) 

where we have used A^L n (— p) = A£*(p), since A M (x, r) is real. After further simplification, 
we get: 

Z = J[dC}[dC}e W {(C,FC)}J[dA»}exp{-^(A»D^A»)}, 

= det(F)(det(ZV)H ; (17) 
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here 



F = i(3p-n, 



= (p 5^ - VvVv + ■ 



;is) 



The determinants are calculated by summing over space-time indices as well as field indices. 
Further manipulation leads to 



\nZ = lndet(/3p-n) — -lndet(D MV ), 
= lndet(/?a> n ) - -lndetfD^), 



with 



r p 2 -pl + ^p -PoPi -P0P2 



D 



-P0P3 



-P0P1 

-P0P2 
-P0P3 



In the temporal gauge, 



det(£> 



p 2 - p\ -pxPi -P1P3 
-P1P2 P 2 ~ p\ -P2P3 
-P\P-i -P2P3 P 2 ~ Pi 



P 2 



a 



(19) 



(20) 



(21) 



Hence 



\nZ 



lndet (^ n )-^lndet(Vp 4 ^) 
2 a 



Tr In {0 2 uo 2 n ) - Tr In (-/?Va;^ 



a 



(iui^+pto+e^^. 



(22) 



n p 



Absorbing the second term in Eq. ( p2| ) into the overall normalization of Z, we are left with 
the expression 

ln ^ = -EEM/? 2 (^ + p 2 )L 



(23) 



n p 



which may be further reduced by applying the formulas |16| : 

r/9V ^2 



In [(27m) 2 + /3V 1 
+00 1 

E 



# 2 + (2vrn) 



+ In [1 + (2vm) 



n 2 + 



2tt' 



2vr 2 / 2 



(24) 



Finally, incorporating the /3-independent term into the normalization of Z, and replacing 
P by 

V 



E 



(2tt) : 



J d 3 p, 



(25) 
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we find that 



In Z = 2V 



IPl 



(26) 



Eq. (ffijp, which is just the usual black-body radiation formula for the two physical gluon 
polarizations, agrees with a similar calculation in the Feynman gauge \TE\. We now proceed 
with the more challenging task of computing the thermodynamic pressure at two loops. 



4 The Q.C.D. Pressure at Two Loops 



Higher-order terms in the partition function may be evaluated by forming bubble diagrams 
and applying the usual Feynman rules. The pure gauge contribution to the pressure at two 
loops is contained in the diagrams of Figure 1. The first of these is called the oyster diagram, 
the second the bowtie diagram. The thermodynamic pressure was the first gauge-invariant, 
physical quantity to be calculated in finite-temperature Q.C.D. at the two loop level. The 
computation was carried out by Kapusta [[17]] in the Feynman gauge who obtained the 
following result: 



pglue 



9 \- v mi 



144' - 9 



where N c is the number of colours, and N g 



N c N g T\ (27) 

N c 2 — 1 the number of gluons. N c and N g 
arise from counting the number of particles participating in the interaction. As stated in 
the Introduction , our aim is to calculate the 2-loop pressure in the noncovariant temporal 
gauge, by using the method of ^-functions, along with a special version of the unified-gauge 
prescription. 

Applying the Feynman rules of Section 2 to the oyster diagram in Figure 1, we obtain: 



poyster 

r 2 



d 3 k 



i i 



(2tt) 3 (2tt) 3 (p + k) 2 p 2 k 2 

S^ip + 2k) T + 5 aT (2p + k)^ - 5 Tfl (p - k) a 

5 a , u (p + 2k) T , - 5 UT '(p - k)a> + <W(2p + k) l 
(p + kf(p + k) a ' 



+ 



(p 4 + h) 2 



TT P T P T 



5 



Pi 
Wk u ~ 
k 2 



(28) 
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where p/4 = 2irnT } k^ = 2irlT } and //, v, r, r', a, a' are only summed over 1,2,3, while n = 
1,2, .. . and I = 1,2, . . .. Multiplication gives rise to a total of 720 terms. Contracting 
indices we find that the resulting expression reduces to 63 tadpole-like integrals, many of 
which are equal to zero. The contribution from the oyster diagram is, therefore, given by 

+oo +oo 63 

P 2 oyster = g*N c N g T £ T £ J> ( 29 ) 

n=— oo l=—oo i=l 

The integrals are summarized in the Appendix. They are related to the familiar 
tadpole integrals which arise in massless theories such as Q.C.D. Let us briefly discuss this 
important class of integrals. 

A typical example, at zero temperature, is the tadpole integral /, 

d A p 1 



w (30) 

In the context of dimensional regularization, such integrals may be formally set to zero 
rj| . The finite-temperature version of (|30| ) reads 



with = 2imT. It is instructive to solve Eq. ( |3TD by using both the contour method and 
the zeta- function method. The contour method gives [T5|, TS]: 



T l77^—f d Po 2 _ 9 -/?coth(^), P=l/kT,k = l, (32) 



or 



(2tt) 3 [ 2m Jioo-e * p 2 - p 2 V 2 erto* - 1 

1 /"ioo+e 1/1 1 

+ / No- 



where 



and 



2iri J-ioo+e p 2 - p 2 v 2 _ 
J VaC +/ matt , (33) 



vac /• d 4 p ( 1 \ . 

1 = / 77T^4 V — ^ — 2"j =0 ' Po = «P4, (34) 
7 (27r) 4 ^p 4 2 + p 2 ' 

matt f d 3 p 1 1 

(27r) 3 ]p[e0lpl - 1' 
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47T f°° p 



/ dp— 3 , (integrating over the angles) 

Jo e^ p — 1 



(2< 

An 1 
rF 2 C 2 - 



(2tt) 3 v ' (3 2 ' 1 
12' 



(35) 



Thus I = T 2 /12. 

Next we compute the tadpole in Eq. (|3~T|) by employing the method of zeta-functions 
|J. Making a Wick rotation to Euclidean space and using the Schwinger representation for 
the denominators, we get 



T Y, J Q dae- ap * 2 J 



2 r d p _.,„j 



n=— oo 



(2tt> 



n=— oo 



^ /n (2tt) 2 ' j - 1 



1 1 + °° 
= 8^ r( ~2 )T S N> ^ = 2 ' ( 36 ) 

n=— oo 

The sum over \p±\ looks divergent, but is actually finite in the spirit of analytic continuation. 
The relevant analytic continuations of the gamma and zeta functions are 



C (l - a) = 7 r- Q 2 1 - Q r(«)C(a) cos (— ), (37) 

T -m+- = i 1 — _v m = l,2, .... 38 

V T l-3-5---(2m-l)' ' ' y 1 

Applying formulas fl37|) and (^) to Eq. fl36|) , we see that 

' = ^TVr(-i) C (-i), 

T 2 

- 12< < 39 > 

since T(-l/2) = -2y/n and ((-1) = -1/12. Note that Eq. © agrees with Eq. ©, and 
that both results vanish in the limit T — > 0. This answer is consistent with the conclusions 
in refs. @. 

There are two other integrals which are similar to tadpoles and arise in the calculation 
of the thermodynamic pressure. They are included here for completeness. The first of these 
integrals is calculated as follows: 



rp g I fV_ g f d 3 P p 4 2 + P 2 



(2tt)3 n ^J (2vr)3p 4 2 + p2 
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n=— oo 
r 3/2 



2T 



7T ' 



(2tt)* 



2 v ^ + 2 v ^F)(27rT) 3 C(-3), 



The second integral may be decomposed as follows: 



where 



while 



Thus 



T E 



d 3 p 

(2tt) s 



P=h + h 



h = T J2 



d 3 p 2 



d 3 p 



(2*) 

oo „ 

W2iy/ (2x)3 

T(2ttT) 2 2C(-2) y "" P 



(2tt) 3 ' 



0, C("2) = 0, 



r E 



n=— oo 
+oo 



d 3 p 2 
(2vr) 3P ' 



-7P 



n=— oo 

+oo 3/2 „oo _ / ,Q ; N 

_ (2^/o ^e-^ 2 ^ 4 2 (-3/2) 7 - 5 / 2 + ^(-5/2) 7 - 7 / 2 



r E 



+oo 



r E 



i r(-3/2) 



^3/2 (p4 2)-5/2 



1 



t E 



d 3 



(27T) 



P 2 n 
p = 0. 



(40) 



(41) 



(42) 



(43) 



(44) 



We now have enough machinery to compute all 63 integrals in Eq. (p9|). Their values 
are listed in the Appendix. Substituting these integrals into Eq. (p9|), we only need to 
complete the indicated summations in order to derive p°^ s * ei '. 
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For instance, 



-oo +00 



t 2 E E 

n=— 00 Z=— 00 



16tt 2 



(n rp\2 00 00 
167T 2 ^ A ^ 

T 4 [C(-1)]" = T 4 [ 



n=l Z=l 

r 1 I 2 

12-1 ' 



The final expression for the oyster diagram reads 



oyster = 2 ^ 
144 



(45) 



(46) 



The contribution from the bowtie diagram in Figure 1 is much easier to calculate. We 
have 



jbowtie 



+00 



-g 2 N c N g T E 



5 xp + 



p x p p 
P 2 , 



+00 



n=— 00 
,2 ^2 



d 3 p „ 

T 

(2tt) 3 


+00 

Z=— 00 " 


r d 3 k 
1 (2tt) 3 


1 

p 2 k 2 










k p k u 
k^ 


), A. 




= 1,2,3, 


d 3 p 

(2tt) 3 


+00 

l=—oo " 


r d 3 k 
1 (2vr) 3 


1 

p 2 k 2 



6 + 2^ + 24 + 24^-2 (p - k) 

P4 



^ 2 ,~ 2 2 ,~ 2 2 / ' 

n-4 ^4 A/4 /^4 A.4 / 



+00 +00 1 ri 

-g 2 N c N 9 T 2 E E t 

n=— 00 i=-oo 



16tt 2 
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(47) 



3 ° ' 9 144' 

Combining Eqs. fl4"6j) and ( f4"?|) subject to their respective symmetry factors, we find that 
the pure-gauge contribution to the two-loop pressure is given by 



pglue _ 1 poyster 1 p bowtie 

r 1 — T7T- r 2 r —-To , 



12 



^-N c N g T\ 



144 



(4£ 



This value agrees with the Feynman-gauge result in Eq. 

5 Conclusion 



In this article we have used the temporal gauge, Aq(x) = 0, to derive the thermodynamic 
pressure in Q.C.D. at one and two loops. Our answer for the one-loop pressure, Eq. (|26|) , 



11 



agrees with the well-known formula for the black-body radiation, while our two-loop result 
in Eq. fl48|) turns out to be identical to the Feynman-gauge answer. 

These results were achieved by replacing the traditional contour method by the method 
of zeta-functions, and by exploiting a special version of the unified-gauge prescription to 
handle the spurious poles of the gluon propagator, Eq. (|2]). These spurious poles may be 
treated by the causal prescription |13| (in Minkowski space): 

1 



Po 
1 



(Po) 2 



^emp =Um £>0; 
Po + ie 

temp =lim (P^_\ 2 e>0 (49) 

which is well defined for po = 2irinT = 0, namely at n = 0. Consistent application of 
this prescription enabled us to evaluate unambiguously all 63 two-loop integrals in the 
Appendix. 
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Appendix 

ovstjGr 

We list here the 63 integrals 19 needed in the computation of P 2 in Eq. 



h = 
h = - 

h = + 

h = + 

h = + 

h = + 

h = - 

h = ~ 

ho = - 
hi = + 

Il2 = ~ 

hs = + 
Iu = + 
hs = - 

/l6 = - 

Ivt = + 
hs = + 



1 



p 4 (P4 + h) h 
Pa 



d 3 p f d 3 k 



(2tt) 3 7 (2tt) 3 ' 



(p 4 + k 4 ) k 4 (k + p) 2 
21 k d 2 



\P4\\h 

(p 4 + k 4 f k 2 (k + p) 2 ^P±l !67r 2 _ 



21^ 



33 p 4 2 



(p 4 + k 4 ) 2 k 2 (k + p) 2 

42 p 4 k 4 
(p 4 + k 4 fk 2 (k + p) 2 
5 k 2 



33 
42 



(g4 ~ 
fc 4 (p 4 — fc 4 



p 4 



) 2 ; 

p 4 - 



(p 4 + k 4 ) p 2 (k + p) 2 

2 



( P4 + h) 2 h 2 (k+p) 2 
>1 



27p 4 2 



27 



■lofi 

k 4 



Pa 



16tt 2 



NN 

16vr 2 

16tt 2 
0. 



(?»4 + fal) P 2 fc 2 

9 

7 2 — = °> 

(p 4 + k 4 ) k 2 

12 k 4 

-^2— = 0- 

p 4 (p 4 + k 4 ) p 2 

(k+p) 2 _ 

p 4 (p 4 + k 4 ) 2 k 4 k 2 
p2 

2p 4 (p 4 + k 4 ) 2 k 4 k 2 
(k + p) 2 _ 
4 (p 4 + k 4 f k 4 2 k 2 " 
(k + p) 2 _ 
2p 4 2 (p 4 + k 4 ) 2 k 2 ~ 

k 2 

4 (p 4 + A; 4 ) 2 A; 4 2 (fc + p) 2 ' 
k 2 

2 = °> 

2p 4 (p 4 + k 4 ) k 4 p 2 



(p 4 + k 4 ) 2 



\p 4 \\k 4 



16tt 2 



0. 
0. 



42 p 4 k 4 



(p 4 + k 4 ) p 2 (k + p) 2 
1 



42 



p 4 (k 4 -p 4 ) 



k 4 



\p 4 \\k 4 



16tt 2 



(p 4 + k 4 ) 2 k 2 



0. 
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( P4 + hY(k + P y 

30 j?4 hi 
(pa + k^f p 2 k 2 

27 k 4 2 _ 
(Pa + A; 4 ) 2 p 2 k 2 

^^ = 0, 

Pa (p 4 + k A ) k 2 

8(k + p) 2 



(Pa + k A ) p 2 k 2 (Pa + k 4 



II 
-IS 
-IS 



p 4 3 k 4 



3 P A 2 (k + p) 2 _ 6 
( P A + kA) 2 k 4 2 p 2 k 2 (pA + k 4 ) 2 k 4 



5 Pa p 2 



(Pa + h) kA k 2 (k + p) 2 



12 pa 3 



(pA + kAYkAk 2 {k + pf 



J 19 = + f f j = °> 

J J Pa 2 (Pa + h) 

j ■ ■ 6 

J 20 — ~~ 

hl = ~ 

122 = — 

-^23 — — 

hi = + 

Its = + 
he 
hi 

hs = + 
hg — — 
ho — — 
hi = - 
I32 = — 
-^33 = _ 
hi = + 
hs = + 
he 
hi 



0, 



-30 



-27 



p A k A 



(Pa + kA) 2 



k& 



(Pa + k 4 ) 



\Pi\\kj 
16tt 2 

\M\k4 
16tt 2 



16- 



p 4 k 4 



\pA\\k 4 



16tt 2 



\Pi\\h 

167T 2 



= -10 



(PA ~ kA)PA 



Pa 



\Pi\\h 
16tt 2 



= 12 



(Pa - k A f 
p 2 k. 



M|fc4 

16tt 2 



7p 4 (k + pf 



14- 



Va 



(pA + kA) k 4 p 2 k 2 (pA + k 4 ) 



II 



12 pa 3 



= -12 



Pa 



4 

\Pi\\kj 
16tt 2 

\pA\\k 4 



(p 4 + kA) 2 k A 



(pA + k 4 ) k 4 p 2 k 2 
6 p 2 = 
(p 4 + kA) 2 kA 2 k 2 
k 2 

4 (p 4 + kAf kA 2 p 2 

(k + pf 2 

4 (p 4 + k A ) 2 k 4 2 p 2 k 2 ~ 3 (p 4 + fc 4 ) 



16tt 2 



0, 



= 0, 



P 4 



|p 4 ||fc4 



167T 2 



P 



22 



4 (p 4 + k 4 ) 2 k 4 2 k 2 (k + p) 2 3 
A; 2 



16tt 2 



4p 4 2 (p 4 + kA) p 2 
,2 



-IS 

-II 



4 (p 4 + kA) 2 k 4 2 k 2 

p2 

^Pa 2 (pa + kA) 2 k 2 



(k + pf 



2 ( P a + A; 4 ) 2 k 4 2 p 2 



0, 
0, 
0, 
0, 



14 



// 



4p 4 2 (Pa + fa) p 2 

.-II 



// 



/ 38 = 

-^39 
ho = 
hi = + 

'« = -JJ 
'« = +// 

'« = "// 



(k + p)' 



P4 2 (P4 + fa) p 2 

12 k/ 



P' 

3-^ = 0- 



Pi (p 4 + fc 4 ) (/C +p) 2 

5 fc 4 A: 2 

Pa (P4 + ^) 2 P 2 (fc + p) 2 
9 
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Figure Caption 

Figure 1. Diagrams contributing to the thermodynamic pressure in Q.C.D 
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